Mean-density Bogoliubov description of inhomogeneous Bose-condensed gases 
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A mean-density description of spatially-inhomogeneous Bose-condensed gases based on Bogoli- 
ubov's method is introduced. The description assumes only a large mean atòmic density and so 
remains vàlid when the mean field collapses due to phase diffusion. A spread in the number of parti- 
des in the condensate is shown to lead to an anomalous coupling between the condensate and excited 
modes. This coupling is due to the dependence of the condensate spatial wavefunction on particle 
number and it could, in principle, be used for reducing particle fluctuations in the condensate. 

PACS number(s): 03.75.Fi, 03.65Bz, 67.90,+z, 05.30. Jp 



Thcre has been a surge of interest in Bosc-Einstein 
condensation following the condensation of alkali gases 
in 1995 H. A defining feature of this recent work is the 
spatial inhomogeneity of the gases due to their confinc- 
ment in relatively small magnètic and optical traps [Q. 
Inhomogeneous Bose condensates also occur in the pres- 
ence of vortices 0j and in the surface region of superfluid 
helium B . Much theoretical work has focussed on mod- 
ifications of Bogoliubov's mean-field method to include 
the inhomogeneity arising, in particular, from quadratic 
trapping potentials. For partides with a repulsive inter- 
action the inhomogeneity leads to two important effccts. 
One is that the mean energy per particle depends on the 
number of partides H|| . Wright et al. [Q and Lewenstein 
et al. [||^] have shown that this can lead to ràpid phase 
diffusion and collapse of the mean atòmic field in a time 
much shorter than the lifetime of the condensate and so 
mean-field descriptions of the condensate in these cases 
remain vàlid only for relatively short times. The other 
effect is that the size of the ground-mode spatial wave- 
function increases with the number of the partides |||(| 
and so changing the number of partides in a condensate 
without corrcspondingly changing the spatial wavefunc- 
tion will produce occupation of excited trap levels. This 
implies that a condensate with a nonzero mean field, and 
thus with a spread in the number of partides, is couplcd 
in some way to the excited levels of the trap. However, 
such particle-fluctuation-dependent couplings are not in- 
cluded in descriptions of condensates to date. The aim 
of this Letter is to introduce a description of a condensed 
gas at zero temperature that takes account of these two 
important effects. The new description depends only on 
the condensate possessing a large mean atòmic density 
without the necessity of a large mean field. Since in other 
respeets it follows Bogoliubov's method it could be called 
a mean-density Bogoliubov description JlCj ] . The signif- 
icant features that set it apart from conventional mean- 
field Bogoliubov descriptions are that it remains vàlid 
even when the mean field collapses and it includes an 
anomalous coupling between the condensate and excited 
modes which depends on the deviation of the condensate 



particle number from its mean value. This coupling could 
be used, in principle, to reduce the particle fluctuations 
in the condensate. The new description has important 
consequences for studies of the quantum-field nature of 
condensates such as state preparation JÏT[ , state recon- 
struction JO], weakly coupled condensates conden- 
sates as open systems 14 and atom làsers |Ï5| . 

We begin with the Hamiltonian for a weakly interact- 
ing boson gas in second quantizcd form 



H = I d 3 r¥(r) \£+ -*t(r)*(r)J vÈ-(r) 



(1) 



where í'(r) is the atòmic field operator at position r, 
C = — t^-V 2 + V(r), V(r) is the trapping potential, 
u = Anh a sc /m, a sc is the s-wave scattering length and 
m is the atòmic mass. The objective is to find an ap- 
proximation to H which is correct to order unity. It is 
convenient to expand , 5(r) in terms of an orthonormal 
basis of spatial modes íf6|JlT 



*W =^S n £„(r) 

n=0 

where à n is the annihilation operator associated with the 
spatial mode £„(r) of the trapped gas with [à n ,àj,J = 
S n ,m, [à n ,à m ] = and J d 3 r£* (r)£ m (r) = 5 n , m . The 
condensate is assumed to be in the "ground" mode rep- 
resented by à n and £ n with n = and so (àjào) — Nq ^> 1 
whereas all "excited" modes are weakly occupied, i.e. 
(àjjà n ) <SÍ Nq for n > 1. Also, the spread in the number 
of partides in the ground mode is assumed to be no larger 
than Poissonian. The effect of the spatial-dependence on 
particle number is most clearly seen by writing the atòmic 
field operator in the form ^(r) = ào£o(r)+ài£i(r)+V>2(r) 
where 

oo 

1p n (r) = ^ «mCm( r ) • 
m—n 

Substituting into Eq. (||) yields 

T~i = 'Hoi + Hrem (2) 
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where 7ioi contains all terms involving óo, àj, ài and à\ 
but not -02 nor 

The most significant contribution to 7ioi are the terms 
involving only the operators ào and à and these are given 
collectively by 



H g = / d 3 ràlC (r) 



C + ^àlàoMr)\ 2 \ àoCo(r) . (3) 



The normalised function £o(r) that minimizes the mean 
of this expression is a solution of the Gross-Pitaevskii 
equation 



[£ + AH£o(r)| 2 Uo(r)=/^o(r) 



(4) 



where N = (àQ^Qàoào)/(ànào) ~ No and n is the Chemi- 
cal potential. Using Eq. (Q) to remove ££o( r ) from Eq. 

(i yicids 



7Y 5 = ^à à 



— -a (a a - 2iV)ao 
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FIG. 1. Spatial mode functions as a function of the radial 
coordinatc r = r|, in units where ro = -jK/rruo = 1, for a 
spherical harmònic trap of frequency 1000Hz with a conden- 
sate containing 10 5 rubidium àtoms (m = 1.44 x 10 - 25 kg, 
a sc = 10 nm). The sòlid and dashed curves are the Lambert 
function S| and Thomas-Fermi |J approximations, respec- 
tively. 



where a n = / <i 3 r ( r ) 1 2rt ■ This result shows that 
the Gross-Pitaevskii equation (Q) defines the minimum- 
energy spatial wavefunction Ço( r ) of the condensate for 
arbitrary states of the condensate (for a given N) in- 
cluding those for which (ào), and thus the mean field 
(*) w (oq)^o, is zero. 

The next most significant contribution to 7i i is given 
by terms which describe a coupling between the n = 
and 1 modes. For example consider the terms containing 
the product £*( r )£o( r ) or its complex conjugate to first 
order, i.e., the terms à\ào J d 3 rÇl(r)£Ç (r) + h.c, which, 
on replacing ££o(r) using Eq. (||), become 

a\à Q f d 3 rfí(r) [/i - NuMr)\ 2 ] £o(r) + h.c. (5) 



The contribution of the first term in the square brackets 
is zero because £o( r ) and £i(r) are orthogonal. If the con- 
tribution of the second term in square brackets was also 
zero then there would be no coupling of this type between 
the ground mode and the excited modes. Therefore the 
part of the function |Ç (i")| 2 £o( r ) which is orthogonal to 
Ço(r) represents the spatial mode function coupled most 
strongly to the ground mode. £i(r) can be made to be 
this maximally-coupled mode by setting 



a(r)=/3[|Co(r)| 2 -a 2 ]e (r) 



(6) 



where (3 = («3 — a|) -1 ^ 2 is a normalisation constant. 
Fig. 1 illustrates the mode functions Ço and £1 for typi- 
cal experimental parametcrs. With this choice for £i(r) 

expression (||) becomes — ^à\ào + h.c. Evaluating the 

integrals in the remaining terms in 7ioi using £i(r) de- 
fined by Eq. (^) and the fact that £o( r ) is a solution of 
Eq. (0) yields 



01 = /ia a + -^-à (aja - 27V)à 



-(a à - N)à\ào + h.c. 



Hià[ài 



+7àíài(2a a — N) + -(àjà{à à + h.c.) (7) 



where 7 = [(3 2 (a4 



72/ 



2a.2\u and fi\ = fi 



|^/d 3 r^[|eo| 2 V 2 eo-V 2 (|eo| 2 eo)]- Terms of order one 
or lower in ào or aj have been neglected because their 
magnitude, being 0(l/^/No), is relatively small. The sec- 
ond term in Eq. (^) produces phase diffusion of the mean 
field [0JÏ7]]. This term is quite distinct from the corre- 
sponding "momentum"-squared term in Lewenstein and 
You's mean- field model B]. In fact, the "momentum"- 
squared term in is essentially a linearized version of 
the above term for the special case where the conden- 
sate has a large mean field and a phase of zero. The 
next term in Eq. (0) represents one of the main results 
of this work: it describes an anomalous coupling between 
the ground mode and mode 1 whose strength depends on 
the particle fluctuations in the ground mode. This term 
embodies the coupling to the excited modes mentioned 
in the opening paragraph. The last three terms represent 
the "usual" terms for excited modes in the sense that re- 
placing ào with gives the corresponding terms in the 
mean field approach. 

We can get some idea of the effect of the coupling be- 
tween the ground mode and mode 1 by considering the 
action of the Hamiltonian term 7ioi alone. This ignores 
the effect of the other excited modes, however, these ef- 
fects should be small if the system is in an approximate 
stationary state of 7i re m in Eq. (^). 7ioi conserves the 
total particle number and so one easy way to follow its 
effect is to consider the evolution of a fixed number of 
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partides, for example, with initially M partides in the 
ground mode and partides in mode 1. The value of 
N need not be equal to M — 1, i.e. £o( r ) need not be 
the lowest energy wavefunction. (The general solution 
for an arbitrary state \A) with a spread in the num- 
ber M of partides in the ground mode can then be ob- 
tained by an appropriate linear superposition of the So- 
lutions for different vàlues of M but the same value of 
77 = (A\àlàlà à \A)/(A\àlào\A).) It is not difficult to 
show for M 3> 1 that the evolution of the mean particle 
number in mode 1 is given to a close approximation by 



(à\ài(t)) = ci sin 2 (u/i) + c 2 [cos(ui't) — 1] 



(8) 



which indicates that partides oscillate between the two 
modes. Here_ci = [(M-y) 2 _ + u 2 (M - N) 2 M//3 2 ]/(huj') 2 , 
c 2 = u 2 {M-N) 2 M[{M-N)(^-im 2 )+ l· i 1 - l· L} 2 /P 2 (huj')' i 
and (Tia/) 2 = [ 1 (2M-N)-(M-N)ua 2 +^i-^] 2 'l 2 M 2 . 
In the Thomas-Fermi regime || the ground mode func- 
tion for a harmònic trap of angular frequency lo is approx- 
imated by £o( r ) = [(/•* — \v\ 2 h,LO /2tq) / Nu] 1 / 2 where r = 
y/h/müj for which /i = BTilo/2, ua 2 — Bhu>2/7N, u/P = 
Bhu2/V67N, j_= BTllo20/77N and ^i = M + 637^/25 
where B = (15Na BC /r ) 2 / 5 . Taking 7 f» ua 2 ~ 2u//3 and 
N ss No leads to the simplifications 

ci « [M + ~(M - iV ) 2 ](a sc /ro) 4 / 5 (287V 1/5 )- 1 
c 2 w (A/-iVo) 2 (16AÍ)" 1 



and lo' ps ^Sò^JM/Nqlü ~ 4cj for r /N a sc < 1- For a 
condensate of rubidium àtoms tightly confincd in a trap 
of frequency 1000Hz the number of àtoms oscillating be- 

4/5 

tween the two modes is of the order of c\ w N /470. 
This can be a significant number for large condensates 
containing Nq ^> 500 àtoms, e.g. Nq = 10 6 gives ap- 
proximately 200 = \/N~n/h àtoms oscillating between the 
modes. For a Poisson-distributed occupation of the con- 
densate a quarter of these àtoms are due to the anoma- 
lous term. 

The dependence of the coefficients C\ and c 2 on the 
difference M — iVo suggests a method for reducing par- 
ticle fluctuations in the ground mode that, despite its 
practical limitations at present, serves to illustrate the 
nature of the anomalous coupling. Imagine that after 
evolving according to Eq. (||) for a time ir/ui' the pop- 
ulation of mode 1 is rapidly depleted by some means. 
Although difficult in practice, in principle it would be 
possible to do this without simultaneously depleting the 
ground mode using the method proposed by Walsworth 
and You |ÏS|] . This removes 4c 2 partides, on average, 
from the system. After repeating this process many times 
the population of the ground mode either approaches A^ 
or decays rapidly to zero depending on whether initially 
M > Nq or M < No, respectively. Evidently this cyclic 
process is a measurement of the occupation of the ground 
mode that distinguishes only between initial occupations 



M that are above and below Nq. Thus, by linearity, if 
the ground mode begins in a pure state which has, e.g., 
a symmetrical spread in particle number M about the 
mean A^o it will eventually be transformed into a state 
which has, with equal probability, either approximately 
A^o or zero partides. This method, therefore, has a 50% 
efficiency in reducing particle fluctuations in the ground 
mode. 

The final task in developing this mean-density Bogoli- 
ubov description is to identify the quasiparticle modes. 
The expansion of the Hamiltonian in Eq. (||) was cho- 
scn to express most clearly the coupling between the 
ground mode and mode 1. It is more appropriate here, 
however, to write the atòmic field operator as í'(r) = 
^o£o ( r ) + "01 ( r ) an d expand the Hamiltonian to second 
order in V>i(r). Transforming off the uninteresting Chem- 
ical potential gives the Hamiltonian in grand canonical 
form, Tif = fi — j d 3 r/i^t (r)^f(r), for which the requircd 
expansion can be written as 



7i — TL' g + Tíge 



(9) 



where H' g , Tí ge and H e are, respectively, independent, 
linear and bilinear in ^i(r) and ^í(r). Choosing £o( r ) to 
minimize the mean of H. g , as before, yields 



H = —a (a a 



2N)à 



(10) 



The foregoing analysis also shows that fí ge represents the 
coupling between the ground mode and mode 1 as given 
by the third term on the right side of Eq. (Q), i.e. 



u 4. — + 
Wge = -g(%a - N)a[a + h.c. 



(11) 



It is useful to consider this coupling term as a perturba- 
tion on a system whose Hamiltonian is given by Tí' g +Tí e . 

In this perturbative approach the quasiparticle modes 
are those that diagonalize 7i e , 



4>l(C - jti)^i + (-■0Mà o ào£o + h - c -) 



2u$l'ip 1 àlà \Ço\ 



Í12) 



The diagonalization can be carried out by generaliz- 
ing the method used by Fetter || to accommodate the 
ground mode operators aj, ào in Eq. (|Ï2| ) which are 
absent from Fetter's mean-field treatment. Since Tíe 
conserves total particle number the required Bogoliubov 
transformation must also conserve particle number. We 
therefore adopt the ansatz |Ï9|| 



L·i(r) 
N 



E 

k=l 



[U k (r)b k 



V k *(r)bl] 



(13) 



3 



where {b k } k are a set of independent quasiparticle opera- 
tors with [bk,bj] = and [&fc,èí] = S k j. Eq. (^3|) implies 
that Í7fc(r) and Vfe(r) are orthogonal to £o( r ); this fact 
can be made more transparent by defining 



U k (r) = u k (r) - ío(r) / d 3 rC (r)u k (r) 



14(r)= Vfe (r)-£o(r) / d ó rC (r)v k (r) 



(14a) 
(14b) 



and working instead with u k (r) and v k (r) which are el- 
ements of the whole space spanned by {£,i( r )}™=o,i,...- 
Since the magnitude of TL e is of order unity one can 
make approximations of relative error of order 
such as replacing the operator à ào with its mean value 
N . Multiplying the first term in Eq. ( |i~2| ) by à^ào/N = 

1 + 0(l/v^Vo), which is effectively unity here, and then 

à \ - 

substituting for —j=ip\ (r) and its hermitian conjugate 

V Ai 

using Eq. ( p^| ) yields an approximate expression for Ti. e 
which does not involve ào nor à . It can be shown using 
an analysis similar to that of Fetter 0] that 



H e — C 



k=l 



fíL·ü k b\b k 



(15) 



to order unity, where uj k , u k {r) and v k (r) satisfy the 
Bogoliubov-de Gennes equations 

Cu k (r) ~ uN\Ç (r)\ 2 v k (r) = huj k u k (r) 
£v k (r) - uN\Ç (r)\ 2 u k (r) = -hu; k v k (r) 



with C = -|-V 2 
C=-EZ=7^ k Jd 3 r[V 3 (r)v* k (r) 



and where 



+ V(r)- l· L + 2uN\Ç (r)\ 2 , 
kj=1 ,^ k fd*r[V 3 (r)v* k (r)+c.c.]/2. 
It is now possible to find the decomposition of mode 1 
in terms of the quasiparticle modes via the overlap 



N 



VN 



Using Eqs. ( p"3| ) and ( |Ï4| ) with the approximate expres- 
sions for u k and v k in the Thomas-Fermi regime found 
by Òhberg et al. [Q for the physical parameters used in 

Fig. 1 gives w 1.755Ò1 - 1.443&J - 0.986(ò 2 - b\), 



N 



that is, mode 1 is composcd of the first two quasiparticle 



modes only |21|. 

In conclusion, this Letter introduces a new description 
of inhomogeneous Bose-condensed gases which is based 
on Bogoliubov's method but assumes a large mean den- 
sity w (àjào)!^! 2 instead of a large mean field. 

The description remains vàlid, therefore, when the mean 
field collapses due to phase diffusion. Eq. @ together 
with Eqs. (10), ( |ïl·| ) and (|Ï5| ) represent its main result: 
an approximate Hamiltonian which is vàlid for arbitrary 
statcs of the condensate and which includes an anomalous 



coupling between the ground and quasiparticle modes. 
This coupling, which has not been described previously, is 
due to the dependence of the ground-mode spatial wave- 
function on particle number. It vanishes on taking the 
mean of 7r with respect to a coherent state of the ground 
mode, i.e. it vanishes in the mean-field approach. The 
mean-density description thus accommodates the two im- 
portant effects of inhomogeneous condensed gases that 
have come to light recently: ràpid phase diffusion and 
the dependence of the ground-mode spatial wavefunction 
on particle number. 

It is a pleasure to acknowledge helpful discussions with 
A. Durrant, S. Hopkins, U. Leonhardt, V. Penna, D. 
Richards and T.B. Smith concerning this work. 
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